A possible new phase of antagonistic nematogens in a disorienting field 
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A simple model is proposed for nematogenic molecules that favor perpendicular orientations as 
well as parallel ones. (Charged rods, for example, show this antagonistic tendency.) When a small 
disorienting field is applied along z, a low density phase JV_ of nematic order parameter S z < 
coexists with a dense biaxial nematic Nf,. (At zero field, A_ becomes isotropic and Nt uniaxial.) 
But at stronger fields, a new phase invariant under n/2 rotations around the field axis, appears 
in between iV_ and Nt,. Prospects for finding the iV+4 phase experimentally are briefly discussed. 

PACS numbers: 64.70.Md, 61.30.Cz, 61.30.Gd 



Practical applications of nematics are based on their 
strong responses to weak external influences. Therefore, 
the phase behavior of nematics in external fields of differ- 
ent symmetries is a subject of many investigations (see 
0,0 and references therein). It was shown that external 
fields can lead to major changes in phase transition be- 
havior. In nematics with positive anisotropy, an external 
magnetic or electric field is 'orienting': it favors positive 
nematic order parameter S. In strong orienting fields the 
first order phase transition between the 'isotropic' phase 
(which for nonzero field acquires finite S > 0) and the 
nematic phase N is suppressed above a critical temper- 
ature. In systems with negative molecular anisotropy, 
similar fields have instead a disorienting effect, favoring 
S < 0. The isotropic phase then acquires uniaxial planar 
orientation; with decreasing temperature or increasing 
density, this phase iV_ has a phase transition into a bi- 
axial phase Nf, (which connects smoothly to N as the field 
is removed). A critical field is expected, above which the 
N—/Nb phase transition becomes second order [[J. 

Previous theoretical investigations of external field ef- 
fects (in cases where the interaction potential was ex- 
plicitly specified 0) were carried out for molecular 
potentials proportional to the second Legendre polyno- 
mial of the relative molecular orientation, i-2(cos W12), 
which is the simplest dependence favouring parallel align- 
ment of nematic particles. In practice, though, real 
molecules have interactions described by more complex 
orientational dependences than this. For example, uniax- 
ial nematogenic particles may display disposition towards 
mutually perpendicular orientations as well the parallel 
one. The molecular interaction then has minima at both 
u>±2 = and u>i2 = tt/2. (Generically these are of un- 
equal depths.) An example of this antagonistic tendency 
is for the case of charged rods. It was shown that the 
Debye-Huckel expression for the electrostatic repulsion 
between two effective line charges at mutual angle LO12 
is proportional to Ismail -1 . As a consequence, two 
charged rods tend to rotate to a perpendicular configura- 
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tion; this interaction is antagonistic to the usual (steric) 
Onsager interaction, favoring parallel alignment, which 
is also present. Such an electrostatic twisting results in 
significant quantitative changes of the nematic phase di- 
agram in zero field (see Ref. Q and references therein). 
Some nonconvex particles, dumbbells for example, should 
have a certain tendency to mutual twisting too. 

The simplest angular dependence of a potential de- 
scribing molecules with a tendency to orient perpendicu- 
larly as well as parallel to each other is the fourth order 
Legendre polynomial 

P 4 (cos W12) = [35 cos(4wi 2 ) + 20 cos(2^i 2 ) + 9]/64 (1) 

In general, this P4 interaction should be taken in addi- 
tion to, not in place of, the customary i-2(coswi2) de- 
pendence. However, absence of a P2 term does not imply 
higher molecular symmetry than when it is present; a 
pure P4 interaction already has unequal well depths for 
perpendicular and parallel interactions. 

The task of this paper is to consider an 'antagonistic' 
nematic in the presence of an external disorienting field. 
The origin of the disorienting field may vary: electric 
or magnetic fields, orienting flow in side-chain polymers 
(with a side group trend to be perpendicular to the back- 
bone), surface influences and so on. The only important 
factor is that the particles want to orient perpendicularly 
to the field direction. Note that besides being ubiqui- 
tous in display applications, which involve both surface 
and electric fields the disorienting effect in a sample 
could be strongly enhanced in some cases by addition of 
impurities [fj. 

Our model comprises a system of uniaxial molecules in 
a disorienting field of strength W, where the potential of 
the molecular interaction with the field is 

v(l) = v{uji) = py(3cos 2 (9i - l)/2, W>0, (2) 

LUi = (8i,(pi) being the orientation of molecule i. The 
above formula assumes that the field is directed along 
the z axis and states that molecular orientations along 
the field direction are energetically unfavorable. The pair 
potential is taken to be a sum of the hard sphere potential 
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for spheres of diameter cr, and of an anisotropic part 

V(l,2)= Wn(-Rl2)Pn(cOSWi 2 ), (3) 



n=1A 



where Wi 2 is the angle between the axes of molecules 1 
and 2; i?i 2 is the distance between their centers; and 



lu n , / , 2 exp(-z„i?i 2 ) 



(4) 



0F/N = p o + J /Hpn/(w) - l]dw 
+ 0W I f(uj)P 2 {cos6)du 



In the limit z„<r — > 0, this model belongs to a family 
of the so-called Kac potentials, for which the mean field 
approach is accurate, and the free energy obeys: 



(5) 



y~] l2f3A n r) / f(uj 1 )P n (cosuJi2)f(uJ2)dujiduj2 
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where r\ = ir(N/V)a 3 /6, f{u>) is a single particle distri- 
bution function, and for the free energy of the reference 
hard sphere system we use the "quasiexact" Carnahan- 
Starling expression, Fq = lnry+(4?7— 3r/ 2 )/(l — 7/) 2 +const. 

To make the mathematics more compact, we give de- 
tails only for the particular case A 2 = 0. (However, 
adding a P 2 (coswi 2 ) term is straightforward.) Using a 
spherical harmonic expansion for P n (coswi 2 ) |j| and not- 
ing that the model under consideration is non-polar, one 
can get the explicit expression for AF = (3F/N — Fo : 

AF = j f{Lj)[]xxf{u)-l]du + pW{P 2 {cos9)) u (6) 

- (3A A ^[(Y 4 ,{u))l + 2 (YM)l + 2 (TZMfj, 

where (•••) u) = J f(u>)(- ■ -)du denotes orientational av- 
erages within a self-consistent one-particle distribution 
function. The superscript 'c' means that the real part 
of the spherical harmonic is taken; © is written as- 
suming that the zx-plane is a symmetry plane. This 
overcomes a continuous degeneracy of the free energy 
connected with the fact that a spontaneous ordering can 
rotate without any energy cost in the xy-plane. 

Minimizing JfjJ with respect to f(u>) we obtain the self- 
consistent single particle distribution function (with C a 
normalization constant) 



f(w) = Cexp{-/W<P 2 (cos0)} 6 
24 



(7) 



+-/M 4 77[<r 40 HLy4oM 
+2 <y 4 c 2 M) w y 4 c 2 M + 2 <r 4 c 4 H) w r 4 c 4 H]}. 



We calculated the phase diagrams of this system using, 
where appropriate, the (standard) common-tangent con- 
struction to find coexisting states of different density. 



The symmetries of all the resulting phases can be denned 
using three order parameters: 



S z 



= (P 2 (cos0)) w , 
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(8) 
(9) 
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Here S z and S x describe alignment of molecules along 
the corresponding axes, whereas S xy describes ordering 
along the y axis as well as along the x axis; when S x = 0, 
S xy describes ordering with fourfold symmetry. 

In Fig. n the phase diagrams are presented for dif- 
ferent values of the ratio W/A4 describing the exter- 
nal field strength in units of the molecular interaction. 
(We set A2 = as explained above.) At zero field 
(W/A4 — > —0) the system undergoes a first order phase 
transition from an isotropic to a uniaxial nematic N; this 
phase has S z > 0, S x = S xy = 0. Small disorienting fields 
(W/Ai = 0.5) transform the system symmetry, so that 
the isotropic phase at low densities is replaced by a uni- 
axial planar phase (this has S z < 0, S x = S xy = 0). 
At intermediate densities (between the bold solid lines 
in the figure) this coexists with a highly ordered biax- 
ial nematic phase Nj,. The biaxial has S z < 0, S x > 0, 
S xy > 0, and becomes a stable single phase at high densi- 
ties. At higher fields the orientational transition becomes 
second order at high temperatures, joining onto the first 
order transition (which is accompanied by the miscibil- 
ity gap) at the tricritical point. The wide miscibility gap 
at low temperatures is connected with the appearance of 
an effective attraction between oriented molecules, which 
causes the gas liquid phase transition. Molecular poten- 
tials of typical thermotropic nematics are characterized 
by strong isotropic attractions, in addition to anisotropic 
interactions. In this case the condensation could arise in 
a non-ordered phase. However, some lyotropic systems 
display a wide miscibility gap between phases of differ- 
ent symmetries. So far, the picture is qualitatively just 
as found previously for a pure P 2 (coswi 2 ) interaction 
At low disorienting field strengths, then, the antagonistic 
tendency of our molecules cannot manifest itself and the 
behavior is that of a conventional nematic system. 

However, at strong enough fields (e.g., W/A4 = 2) the 
situation differs strongly from the conventional case: two 
consecutive orientational transitions now take place. The 
first one is always of second order and meets the binodal 
at a critical end point. The resulting phase is invariant 
under rotations by 7r/2 around the field direction; it has 
S z < 0, S x — 0, and S xy > 0. In accordance with the 
general classification this new phase should be de- 
noted as iV +4 . The second orientational phase transition 
is from A^ +4 into N},; this changes from first to second 
order at the tricritical point as before. The new A^ +4 
phase is (like iV&) a result of spontaneous ordering, and 
is characterized by long-range correlations of Goldstone 
type. Note that the molecular potential under investi- 
gation does not possess the symmetry of the A^ +4 phase 
(even for the particular case studied, of A 2 = 0). In 
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FIG. 1: The phase diagram of the system (A2 = 0) in (a) 
the temperature-density (T* = ksT/Ai, r\ = ivpa 3 /6) and 
(b) pressure-temperature coordinates (P* = (iter 3 / &)(P / A4)) 
for different values of the disorienting field (W/A4). N+4 is 
the phase with S xy 7^ 0, S x — 0, Nt denotes the biaxial 
nematic phase. Dashed lines correspond to the the second 
order phase transition into jV+4 phase, thin solid lines likewise 
into Nb; thick lines denote the first order transition. Circles 
are tricritical points, diamonds denote critical end points. 



some respects the situation is similar to that found by 
computer simulations, when a cubatic state appears in a 
system of hard truncated spheres 0] . The iV+4 phase is 
not cubatic but might be called 'square-atic', since it has 
fourfold rotational order in the transverse plane only. 

From the calculated phase diagrams (especially in T— r\ 
coordinates) one can see a nonmonotonic efect of the ex- 



TABLE I: Coordinates of the tricritical point (A% = 0) at 
different values W/A4 of the disorienting field. (Notation as 
inFig.[T]) 



W/Aa 


T* 


V 


p* 


1 


2.02 


0.62 


39.1 


2 


1.70 


0.48 


8.59 


5 


1.30 


0.32 


1.52 


50 


1.52 


0.29 


1.40 



ternal field on the -/V +4 region. The widest ^+4 region 
arises at intermediate field values. We attribute this fact 
to the disorienting field's ability to restrict the part of 
orientational phase space in which particles mainly re- 
side. Indeed, in the limit of an infinite field, the particle 
orientations are confined in the plane perpendicular to 
the field. In this case the system can be described by 
the plane rotator model this has the same phase 
diagram (independent of n) for all pair potentials pro- 
portional to cos(nui2). In contrast, with an unrestricted 
orientational phase space (that is, at zero field), the tem- 
perature at which spontaneous ordering occurs decreases 
with increasing n. Disorienting fields (pushing the sys- 
tem towards the plane rotator state) promote the order- 
ing with S xy 7^ more strongly than that with S x ^ 0, 
so that the -/V+4 phase grows leaving behind competing 
ordered phases such as N b . However, in strong enough 
fields the susceptibility of the order parameter S xy to the 
field decreases, whereas the Nb ordering is not saturated. 
As a result, in strong enough fields the Nb phase begins 
to displace the N + 4 region. 

This field effect on the N + 4 region is closely linked with 
the tricritical behavior (Table[IJ . The nonmonotonic field 
effect on where the tricritical point lies can be explained 
by the existence of two competing tendencies |9|. First, 
the external field furthers the orientational ordering of 
particles and, therefore, causes a more effective attraction 
between particles. This favors the phase separation and 
raises the 'binodal' in temperature. The second tendency 
takes place if the susceptibility of the low-density phase 
(AL|_4) to the field is larger than that of the coexisting 
dense (Nb) phase. This decreases the energy gain of the 
phase separation. The second tendency is very strong 
for our system. Only at very strong fields (W/A4 = 50) 
do the susceptibilities of the coexisting phases become 
comparable, allowing the first tendency to win at last. 

To study the robustness of the predicted A^ + 4 phase 
with respect to changes of the intermolecular potential 
we calculated phase diagrams for A2/A4 = 0.1; 0.2; 0.3 
at W/A4 = 5; 10; 50; 00 in a density interval sensible 
for liquids (r] < 0.6). In Fig. |3 the phase diagram for 
A 2 = OAA4 is presented (setting W/A4 = 5). With this 
potential, the parallel orientation of molecules is more en- 
ergetically favorable than in the case A2 = 0. As a conse- 
quence, the Nb transition temperature increases, whereas 
the temperature of the A^ + 4 transition is not affected sig- 
nificantly. The N + 4 region is reduced in size, but the 
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FIG. 2: The phase diagram of the model J3J with A 2 = 0.1 At 
in the temperature-density and pressure-temperature coordi- 
nates for a moderate disorienting field (W/A4 — 5). (Notation 
as in Fig. 0) For comparison, the location of the iV+4 phase 
for the case A2 = is shaded. 

phase remains stable. The phase transition N+±/Nb is 
first order now (where previously it was second order, for 
this value of W/A4) which should lead to a jump in the 
elastic constants, for example. For all the investigated 
cases with A 2 > the phase transition into Nb is first or- 
der (even at infinite field) whereas the N_/N + 4 transition 
remains second order. It follows from our results that the 
N +4 region is broader in density and temperature at the 
moderate field W/A4, = 10. Moreover, for A 2 = 0.2A 4 
(when the parallel orientation of molecules is 4.4 times 
more energetically favorable than the perpendicular one) 
the iV+4 phase exists at W/A4 — 5; 10; 50, but is absent 
in infinite and zero fields. When the parallel orientation 
is almost 6 times more favorable than the perpendicular 
one (A 2 = O.3A4) the Nb phase finally excludes N +4 . 

In conclusion, we have studied a model system of uni- 
axial particles with antagonistic molecular interactions in 
the presence of a disorienting field causing an easy-plane 
anisotropy. This may lead to spontaneous orientational 
ordering of two types. The first is a biaxial nematic Nb 
with a nonzero value of the conventional biaxial order 
parameter S x — (sin 2 6 cos(2^p))^ (as well as of the uni- 
axial order parameter S z = (P2(cos#)) w ). The biaxial 
order parameters is however zero in the second sponta- 
neously ordered phase, named iV+4, which we predict to 



arise at intermediate density. The symmetry of this phase 
is defined instead by a nonvanishing higher orientational 
average S xy = (sin 2 0cos(4<p)) . This corresponds to a 
cross-shaped distribution of preferred molecular orienta- 
tions, which can be rotated in the xy-plane without any 
energy cost (just like the usual transverse director in the 
Nb phase). This feature is connected with Goldstone- 
mode fluctuations of particle orientations in the plane. 

Device applications for nematic phases are of 
widespread importance p|, and a number have already 
been envisaged for the biaxial nematic phase 0|. But 
thermotropic biaxial nematics have proved to be very 
elusive. One of the causes is that they are displaced by 
other structures. It makes it necessary to study all phases 
nearby (like N +4 ) and conditions of their existence. An 
important question that arises is how to look for the new 
N +4 phase experimentally. N +4 could be confused with 
Nb in its elastic behavior. Both should demonstrate elas- 
tic resistance to the following distortions: 1) splay and 
bend in the xy-plane, 2) twist in the z direction. In con- 
trast to Nb and iV+4, the uniaxial phase A_ does not 
posess such elastic behavior. The new symmetry of the 
N +4 phase also should influence the behavior of defects, 
as is also true for biaxial nematics ^3|. O n the other 
hand, since it has a fourfold symmetry axis, the N +4 
phase is uniaxial with respect to second rank tensors, 
and therefore not distinguishable by simple dielectric or 
optical measurements from an ordinary uniaxial nematic. 

Finally it is worthwile to note that the N +4 symme- 
try could in principle appear without external fields, in 
molecules of a more complex geometry. From this point 
of view the investigation of systems of biaxial particles 
in search of the iV + 4 phase represents an interesting 
prospect. The ordering of long molecular axes eventu- 
ally leads to the restriction of orientational phase space of 
the molecular short axes without any external field. This 
can influence their further orientational ordering with the 
same consequences as discussed above. 

R. O. S. thanks NATO and the Royal Society for fi- 
nancial support. 
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